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ABSTRACT 

We report 350 and 230 GHz observations of molecular gas and dust in the starburst nucleus of NGC 
253 at 20-40 pc (r'-2'') resolution. The data contain CO(3-2), HCN(4-3), C0(2-l), 1^00(2-1), 
C^^0(2-l), and continuum at 0.87 mm and 1.3 mm toward the central kiloparsec. The C0(2-l) 
size of the galaxy's central molecular zone (CMZ) is measured to be about 300 pc x 100 pc at the 
half maximum of intensity. Five clumps of dense and warm gas stand out in the CMZ at arcsecond 
resolution, and they are associated with compact radio sources due to recent massive star formation. 
They contribute one third of the CO emission in the central 300 pc and have ^^CO peak brightness 
temperatures around 50 K, molecular gas column densities on the order of 10^ Mq pc~^, gas masses 
on the order of 10^ Mq in the size scale of 20 pc, volume- averaged gas densities of nu^ ~ 4000 
cm~^, and high HCN-to-CO ratios suggestive of higher fractions of dense gas than in the surrounding 
environment. It is suggested that these are natal molecular cloud complexes of massive star formation. 
The CMZ of NGC 253 is also compared with that of our Galaxy in CO (2-1) at the same 20 pc 
resolution. Their overall gas distributions are strikingly similar. The five molecular cloud complexes 
appear to be akin to such molecular complexes as Sgr A, Sgr B2, Sgr C, and the / = 1.3° cloud in 
the Galactic center. On the other hand, the starburst CMZ in NGC 253 has higher temperatures and 
higher surface (and presumably volume) densities than its non-starburst cousin. 

Subject headings: galaxies: starburst — galaxies: ISM — galaxies: individual (NGC 253) 



1. INTRODUCTION 

NGC 253 is among the nearest and infrared-brightest 
starburst galaxies in the sky and is therefore one of the 
best targets in the local universe for starburst studies. Its 
active star formation is located in the central half kpc of 
this barred spiral galaxy at the distance of 3.5 Mpc {1" 
= 17 pc; Rekola et al. 2005). The starburst luminosity 
of -3 X lO^^Lo (iTelesco HarpedlToSOh suggests a star 
formation rate of ~5 Mq yr~^ according_to_the conver- 
sion from LiR to star formation rate bv lKennicuttI (|^98). 
The intensity of this starburst is approximately 20 Mq 
yr~^ kpc~^ in star formation surface density, among the 
highest (and comparable to that in M 82) in th e sample 
of nearby circumnuclear starbursts in iKennicutt. (^1998i) 
when mergers are excluded. 

Previous high-resolution observations of star form- 
ing interstellar medium in the nucleus of NGC 253 
showed that a quantity of molecular gas, about a few 
1O^M0, in the starbur st region is w arm, dense, and 
kinem a tically disturbe d ([Canzian et al. 1988; Carlstrom 
19901: iPa glione et al." '1995', '2004'; Peng et al. 1996; 
Sorai et a l. 2000; Garcia-Burillo et al. 2000; Ott et al. 
2005; Sakamoto et al. 2006, hereafter S06; Minh et al. 
2007; Knudsen et al. 2007). However, the spatial resolu- 
tions of these observations, including our own (|S06f ), were 
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modest (^3^^) and their frequency coverages limited (< 
230 GHz). 

In this paper we report f -2'' (20-40 pc) resolution 
imaging of the starburst nucleus of NGC 253 in five 
molecular lines and continuum emission in the 230 GHz 
and 350 GHz bands. Several molecular cloud complexes 
are detected in the central 300 pc where the starburst is 
most intense. The gas clumps are associated with known 
signs of massive star formation in the region, suggesting 
them to be the natal clouds of, or a formative environ- 
ment for, massive clusters. The goals of this paper are 
to provide an updated overview of the starbursting cen- 
tral molecular zone and to characterize the star-forming 
molecular cloud complexes there on the basis of our new 
high-resolution multi-wavelength data. The high resolu- 
tion enabled us to directly constrain gas properties, e.g., 
temperature from the line brightness temperatures and 
volume- averaged gas density from the mass and size of 
the molecular clumps. We also compare the center of 
NGC 253 with the center of Milky Way using C0(2-l) 
data of the same linear resolution. We point out morpho- 
logical similarities as well as differences in the gas prop- 
erties between these starburst and non-starburst galactic 
centers. 

2. SMA OBSERVATIONS AND DATA REDUCTION 

The observation s were m ade with the Submillimeter 
Array (SMA; H o et al.ll2004l FI as logged in Table E We 
simultaneously observed three CO isotopologues, CO, 
^^CO, and C^^O, in their J=2-l transitions in the 

^ The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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230 GHz band. We collected 5.7 hr of data toward 
a=00^47"^33n7, (^=-25°iri7Vl (J2000) at the baseline 
lengths between 8 and 509 m f rom six nights of obser- 
vations. Our old data used in lS06l were not combined 
to these data in order to verify our previous results 
with an independent dataset. The 345 GHz observa- 
tions were made in two nights in an array configuration 
called compact-north covering baseline lengths of 7-124 
m. The CO (3-2) and HCN(4-3) lines were observed for 
one night each. Three positions, the center given above 
and two positions 14'' from it along the position angle of 
60°, were observed. The resulting mosaic covers almost 
the same area as the 230 GHz data. Both bands have a 
common flux scale tied to the primary calibrator Uranus 
(IGriffin fc: OrtonI [19931) . We used quasars J0132-169 
and J2258— 279 for our gain calibration and such bright 
sources as Saturn, Uranus, J1924— 292, and 3C454.3 for 
our passband calibration. Pointing was checked during 
the tracks using J1924-292. 

The data were reduced in the same manner as 
in lS06l — initial calibration with the SMA version 
of MMA ([Scoville et ap Il993f ). imaging in MIRIAD 
( Sault et alJ [1995[), and data analysis in NRAO AIPS 
( Bridle Grd^[l99l . Special care was taken to en- 
sure consistency in flux scale between tracks, by compar- 
ing line and continuum intensities on common baselines. 
The uncertainty of fiux calibration is 10% and 15% for 
230 GHz and 345 GHz, respectively. The line data were 
continuum-subtracted in the u-v domain by using contin- 
uum made from the channels off the CO and HCN lines. 
The only other line feature that we noticed in our pass- 
band was a weak HNCO(10o,io-9o,9) line (/rest = 219.798 
GHz) between 1^00(2-1) and C^^0(2-l). The line was 
detected at 5-10 a in a few of channel maps but its total 
fiux is 6% of the C^^ 0(2-1) fiux and its inclusion to our 
continuum affects our continuum fiux by only 2%. Con- 
tinuum maps were made by averaging both sidebands; 
their mean frequencies (wavelengths) are 225 GHz (1.33 
mm) and 345 GHz (0.87 mm). All maps were corrected 
for the sensitivity pattern of the primary beam and mo- 
saicking. Velocities in this paper are with respect to the 
local standard of rest and are defined in the radio con- 
vention. 

The parameters of our data sets are summarized in 
Table [21 Among them are fractions of single-dish fiux 
as detected with our inter ferometric observations. We 
detect most fiux in the galactic center with slightly better 
recovery in the 230 GHz band. Fig. [T] shows the spectra 
and total fiuxes of the lines detected at the center of NGC 
253. 

3. RESULTS AND DISCUSSION 

3.1. Overview of the Central Molecular Zone 

Our line and continuum maps are shown in Figures [21 
[31 [H and [5l The data clearly show the bar-like region of 
about 30'^ (500 pc) length w here the emission has been 
known to be strongest (e.g., S06). Following the con- 
vention for our Galaxy, we call the region the central 
molecular zone (CMZ) of NGC 253. In the core of the 
CMZ are several peaks or clumps of molecular line and 
continuum emission. In particular, there are five strong 
peaks in the 1.3 mm continuum map as listed in Table 
[31 Each peak has a corresponding feature in our 0.87 



mm map, in our CO (2-1) map, and in many cases 
in other molecular line maps presented here. The five 
peaks also have counterparts in existing observations at 
other wavelengths. Some of these clumps have been seen 
in previous interferometric studies of molecular gas at 
lower resolutions (they are easier to see in more recent 
observations and discernible at least in Knuds en et al. 
20071 : iMinh et al."2007': ISa kamoto et al.h2006: .Ott et al 
20051: [Paglione et al. 2004; IGarcia-Burillo et all 120001 : 
Peng et al. 1996). We put our main focus on these con- 
tinuum and line peaks and describe their properties in 
the following. The millimeter peaks are hereafter referred 
to by their numbers in Fig. [2fc) and Table [3l 

3.2. Continuum Properties 
3.2.1. Emission Mechanism inferred from Spectral Index 

The majority of continuum emission at around 1 mm is 
most likely thermal dust emission, judging from the spec- 
tral index between 0.87 and 1.3 mm. The fiux density 
ratio of the 0.87 and 1.3 mm continua is approximately 
3 according to a data comparison with matching range 
of u-v lengths. (The ratio i^o. 87/1.3 is in the range of 
3.2 ±0.5. See Table [3l) This steep increase of fiux density 
toward shorter wavelengths, which was also seen between 
3 and 1.3 mm (S06 and references therein), is consistent 
with the dominance of thermal dust emission at A < 1 
mm. The minimum fraction of dust emission is 0.8 at 
0.87 mm for the fiux density ratio of i^o. 87/1.3 = 3, dust 
emissivity index of /3 = 2 (|Dunne fc Eales|[2001 f). and a 
fiat spectrum contribution from optically thin free-free 
emission. The ratio at the nucleus, i.e., the peak No. 3, 
is i?o.87/i.3 = 3.0 and is not distinctly different from the 
ratio at the other peaks. 

3.2.2. Mass and Column Density from Continuum 

The peak column density of gas at the continuum peaks 
is estimated to be (1-4) xlO^ Mq pc~^ from the 1.3 mm 
emission and to be (0.2-2) xlO"^ M© pc-^ from 0.87 
mm (Table [3l). These estimates use the mass opacity 
coefficient of Hildebrand (1983), /3 = 2, gas-to-dust mass 
ratio of 100, and a dust temperature of 50 K. The peak 
column densities from 1.3 mm data are larger by a factor 
of ~2 than those from 0.87 mm data partly because of 
the larger beam size (and dilution) of the latter data. 
Other possible reasons for any discrepancy between the 
parameters derived from 0.87 mm and 1.3 mm data are 
deviation of dust temperature and j3 from the adopted 
values, different u-v coverages, small contamination to 
the continuum from non-thermal or line emission, and 
our ffux calibration errors. 

The overall uncertainties of our estimates for the peak 
surface densities and masses are about a factor of 5, 
considering the uncertainties of the opacity coefficient, 
gas-to-dust ratio, dust temperature, and calibration of 
our observations. The l argest source of erro r is in the 
dust opacity coefficient. iPollack et al.l (|1994l ) estimated 
its very conservative error to be a factor of 4 for a va- 
riety of grain properties. The gas-to-dust ratio that 
we ado pted is from the rat io in the solar neighborhood 
(^150: [S"odroski et al.lll997h corrected for the metallicity 
at the center of NGC 253, which is 0.1-0 . 2 dex higher 
than in our vicinity (jPilvudn et al.l 12004 l2006 f). The 
dust temperature in the calculation above is from the 
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peak brightness temperature of the ^^CO lines at the con- 
tinuum sources. It is also consistent, within a factor of 
2, with the dust temperatures inferred for the nuclear re- 
gion by fitting the spectral ene rgy distribu tion (SEP) in 
the fa r-infrared (Rad ovich et aLl l2QQl: Me lo et alJl2QQ2l : 
IWeifi et al. 2008). We did not subtract non-dust emis- 
sion from the continuum or correct the continuum flux 
density for missing flux because the two corrections are 
less than a few 10% each and are in opposite directions. 

3.2.3. Association with Known Sources 

The continuum peaks at A ^ 1 mm are spatially well 
correlated with those in centimeter wavelengths, as il- 
lustrated in Fig. [2] (c). Four of the five peaks of 
dust emission have one or more counterparts of centime- 
ter compact sources, in particular the bright centimeter 
pea ks in the nuc lear region customarily called THn af- 
ter [Turner fc H d (1985). The brightest 1.3 mm peak. 
No. 3 in Fig. [2] (c), is associated with the brightest 
centimeter source TH2, which is often presumed to be 
the radio nu cleus and is within 2^^ of the stellar kinemat- 
ical center (|Miiller- Sanchez et al.|[2QTQl ). This millime- 
ter source is elongated along the position angle of ~60° 
in agreement with the elonga ted distributi o n of com- 
pact radio sources obs e rved b y [Turner ([T985h and 
lUlvestad fc Antonuccil (|1997[ ) in the immediate vicinity 
of TH2. Since the radio sources comprise supernova rem- 
nants and H II regions, as well as the possibly active 
nucleus of low luminosity, the coincidence of centimeter 
and dust emission peaks suggests that the 1 mm peaks 
are dust (and gas) condensations heated by those sources 
(i.e., mostly massive stars). 

In addition to the centimeter- wave conti nuum sources, 
two II 2O masers in the galactic center (iHenkel et al.l 
l2004t iHofner et al.l 2006) are spatially associated 
with two of our continuum peaks. Also, the most 
luminous infrared source in the ga lactic center, at- 
tributed to a super st ar cluster (|Keto et al.l 119991 : 
iKornei fc McCraHvl 120091 ), is at the millimeter emis- 
sion peak No. 4 according t o its spatial association 
with the radio source TII7 CKalas fc Wvnn-Williamsl 

[1994: iDudlev fc Wvnn-Williamsl [1991 

iFernandez-Ontiveros et al.l I2009D . Our 1.3 and 0.87 
mm continuum sources are also associated with peaks in 
the molecular line emission described below. 

3.3. Line Emission Properties 

Molecular line maps in Figs. [3] and H] show peaks cor- 
responding to the five dust continuum peaks, less promi- 
nent gas clumps, and extended line emission. The ob- 
served parameters of the line emission peaks are listed in 
Table 11 

3.3.1. CO Excitation and Gas Temperature 

We first note that the peak integrated intensities as 
well as the peak brightness temperatures are approxi- 
mately the same in the two transitions of ^^CO at the 
five most prominent peaks of dust and gas emission. This 
suggests that the bulk of the ^^CO emission from the 
line-emission peaks, or molecular cloud complexes, is op- 
tically thick and is nearly thermalized at least to the 
J = 3 level. Under this condition of (near) local thermal 
equilibrium (LTE), the brightness temperature of an op- 
tically thick line gives the gas temperature at the cloud 



photosphere multiplied by the beam filling factor of the 
emission. 

The peak brightness temperature (Tb) of the ^^CO 
lines, shown in Fig. 21 has maxima of about 50 K at or 
in the vicinity of the integrated intensity peaks. These 
beam-averaged brightness temperatures set a lower limit 
to the beam-averaged temperature of molecular gas at 
each location. Our data therefore suggest that molec- 
ular gas at and around four of the five 1 mm peaks 
has a minimum beam- averaged temperature of 50 K. 
This temperature is not much higher than our previous 
lower limit of 40 K based on the ~3'' resolution obser- 
vations of S06 , suggesting that the extent of individual 
CO-emitting cloud complexes is about l''-3'' (20-50 pc). 
An exception to the ~50 kelvin Tb among the five peaks 
is the galactic nucleus, i.e., peak No. 3, where peak Tb is 
about 30 K in both ^^CO lines. This is quite likely due 
mainly to the fast gas motion around, and large veloc- 
ity gradient across, the dynamical center since gas has a 
small filling factor and a large beam dilution in such a 
location. It is also notable that the region with peak Tb 
higher than 30 K, which is about 3 times higher than the 
peak brightness temperature s of giant m olecular clouds 
in the disk of our Galaxy (iScoville et al.lll987[ ), has an 
extent of 300 pc by 50 pc {17" x 3'') in Fig. |l(a). The 
high temperature in the CMZ of NGC 253 is therefore 
not limited to the most prominent emission peaks. This 
notion is strengthened when we consider that the beam 
filling factor of gas is likely lower at locations further 
away from the peaks of integrated intensities. 

3.3.2. Gas Golumn Density, Mass, and CO Opacity 

The peak column density of gas at the molecular cloud 
complexes is estimated to be ~10^ Mq pc~^ from the line 
emission using multiple methods that we describe below. 
This is in agreement with the estimate from dust contin- 
uum in ^3.2.21 The gas mass of each molecular complex 
is therefore on the order of 10^ Mq for the clumps iden- 
tified in our ^20 pc resolution data. 

Our three methods to estimate gas surface density 
from the line data are the following. Firstly, the sur- 
face density of molecular gas is Smoi = 1 x 10^ Mq pc~^ 
for the ^^CO integrated intensity /i2co = 5 x 10^ K 
km s~^, a typical peak value at the five peaks in both 
CO(3-2) and (2-1), for the ^^CO-Hs conversion factor 
of Xco = 1 X 10^° cm-2 (K km s"^)"^ and the factor 
of 1.4 correction for heavy elements. Here we did not 
make a distinction between the ^^CO lines that we ob- 
served and CO (1-0) for which Xco is most often esti- 
mated, since E'^co{i-q) should be approximately equal 
to /i2co(2-i) and T2co(3-2) foi" the excitation condi- 
tion that we suggested most likely. The Xco factor is 
usually estimated to be in the range of (0.3-3) x 10^^ 
cm~^ (K km s"-*^)"-^ for large metal rich galaxies (e.g., 
lArimoto et al.l 1 1996 1 and references therein) . For the 
center of NGC 253, Xco has been estimated to be 4, 
1.7 1.3, and . 3 in u nits of 10^" cm"^ ( K km s'M'^ 
bv llsrael et all (I1995D. TWeifi et al.l (l2008f) . iMartm et aP 
(2010), and Mauersb erger et al.l (|1996[ ). respectively. 
Secondly, the peak gas column densities at the molec- 
ular clumps are estimated to be Smoi ^ 10"^-^ pc~^ 
from /i3co(2-i) ^ 800 K km s"^ and Smoi ^ 10^-^ 
Mq pc~^ from /ci80(2-i) ^ 300 K km s~^ using the 
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conversion relations from line flux to gas column den- 
sity Smoi = 10^-^/i3co(2-i) and Smoi = 10^-'^/ci8O(2-i), 
where Smoi is in Mq pc~^ and the line flux K km s~^. 
Assumptions used here are that the lines are optically 
thin and from gas in LTE at 50 K and that the molecu- 
lar abundances are [l2CO/H2]=10-^ [^^C/^^q^SO, and 
ji6q^18q]^3qq. ^Yie latter two are from t he lower lim- 
i ts in the center of NGC 253 obtained by iMartm et al.l 
(|2Q1Q[ ). Thirdly, the peak gas column density for 
^ci80(2-i) ^ 300 K km s-i is S^oi ^ 10^- ^ Mq pc"^ if 
we use the relation Hmoi = 10^"^^ci8O(2-i) (Stu tzkilflQQOl : 
iMauersberger et al.llTMl iWild et a l. 1992). Here the as- 
sumptions are that the C^^O line is optically thin, that 
the gas is warm and moderately dense (i.e., Tkin > 40 
K and 10^ < < 10^"^ cm~^), and that the molecu- 
lar abundances are the same as the one assumed above. 
In this third method, LTE is not assumed. Indeed, the 
factor of 2 lower estimate than the LTE value is because 
this method takes into account the difficulty of populat- 
ing levels higher than J = 2. 

The optical depths of CO lines are 15, 0.17, 0.05, and 
25 for i2C0(2-l), 1^00(2-1), 0^^0(2-1), and ^^C0{3- 
2), respectively, for a 50 K gas in LTE with a line width 
of 100 km s~^ and the surface density of Smoi = 

1 X 10^ 

Mq pc~^ (i.e., H2 column density Nb_2 = 5 x 10^^ cm~^). 
These calculated opacities should be representative of the 
line opacities of the cloud complexes because the param- 
eters used are the ones we estimated or adopted for the 
cloud complexes. The above opacities for ^^CO and C^^O 
are indeed comparable to the observationally estimated 
mean opacities of these lines in Table |H Unless ^^C and 
^^O are significantly overabundant than we assumed, the 
^^CO lines that we observed in the molecular cloud com- 
plexes are optically thick. The 1^00(2-1) and 0^^0(2- 
1) lines are optically thin on average although there may 
well be small regions of high column density where these 
lines are opaque, because the gas distribution in each 
molecular complex is almost certainly non-uniform. 

3.3.3. Line Profiles and Self- absorption of CO 

The line profiles of the five most prominent molecu- 
lar complexes are shown in Fig. [6l The line widths are 
around 100-200 km s~^. The line shapes differ from one 
complex to another and sometimes from one line to an- 
other for the same molecular complex. The broad line- 
widths and the variation of line shape among the cloud 
complexes must be partly due to the galactic rotation and 
the large-scale non-circular motion in the barred galaxy 
as seen in the mean velocity map (Fig. [5]) rather than 
due solely to velocity dispersion of molecular gas. 

Regarding the variation of line profile between different 
lines, we note the cases of the molecular complexes No. 
2 and 3. In the former the ^^C0(2-l) and (3-2) profiles 
have a less sharp peak than the ^^CO and C^^O profiles. 
A dip is in ^^CO profiles (or at least in the J=2-l pro- 
file) at the peak velocity of ^^CO and C^^O. In the latter 

CO (2-1) has a sharp dip that is not seen in other lines. 
We suggest these to be self- absorption of ^^CO, since, as 
we evaluated above, ^^CO is optically thick and ^^CO 
and C^^O optically thin. Since the continuum emis- 
sion has peak brightness temperatures of less than 2 K 
at these locations, any absorption much deeper than 2 



K must be mostly toward brighter CO emission rather 
than continuum. The self-absorption then means that 
there is colder gas in front of warmer gas along our line 
of sight. Each absorber can be either inside or outside 
of the cloud complex but is probably within the central 
molecular zone considering the 78° inclination of NGC 
253. The waker (or lack of) absorption in CO (3-2) can 
be partly because the foreground gas is not warm and 
dense enough to have a significant population at J=2. 

An important implication of the self-absorption is that 
the peak brightness temperature, column density, and 
mass of the molecular cloud complexes will be underesti- 
mated using ^^CO for the absorbed gas. For example, the 
molecular cloud complex associated with the millimeter 
continuum peak No. 2 may have the largest column den- 
sity among the five complexes. It is the brightest among 
the five in the optically thinnest C^^O line, even though 
it has comparable or less integrated intensities in ^^CO. 

In addition to the self-absorption, there are probable 
cases of cloud overlapping along our line of sight. For 
example, the molecular complexes No. 4 and 5 have 
double-peaked line profiles. These are likely due to mul- 
tiple velocity components because the double peaks are 
visible even in ^^CO and, in the case of the complex No. 
4, C^^O. The CMZ must be a thin disk-like volume even 
though the most prominent dust and gas features are on 
an one-dimensional string on the sky from our near edge- 
on viewpoint. Thus there must be substantial overlap of 
molecular clouds along our lines of sight. 

3.3.4. Gas Density and Line Ratios 

The volume- averaged gas density is estimated to be 
^ 10^-^ cm~^ for a 20 pc diameter region around 
the molecular emission peaks. This is from the peak 
gas column density of 10^ Mq pc~^ and an assumption 
that half of the column density is from within the 20 pc. 
This density is almost as high as the critical density for 
CO excitation for the J = 2 level, 10"^ cm"^ for 50 K. 
Since ^^CO excitation has additional help from photon 
trapping, it is almost certain that ^^CO is thermalized 
at J ~ 2 as we inferred from the CO(2-l)/(3-2) ratio. 
Our volume- averaged gas density i s comparable to the 
density of n^i^ - 10^-^ cm-^ that IGiisten et all (|2006[ ) 
estimated for the CO-emitting molecular gas in the cen- 
tral 15'' (250 pc) from LVG analysis of multi-transition 
CO data. This agreement between macroscopic and mi- 
croscopic gas densities suggests that the CO (2-1) emit- 
ting gas has a large volume filling factor in the molecular 
cloud complexes. When the molecular gas containing 
^^CO almost fills the volume of each molecular emission 
peak, the effect of beam-dilution must be small, unless 
there is a large systematic velocity gradient at the emis- 
sion peak. Thus the peak brightness temperatures of ~50 
K that we measured in ^3.3.11 must be a good estimate 
of the mean gas temperature at the ^^CO photosphere 
(with the likely exception of the molecular peak No. 3 
at the galactic dynamical center.) 

The intensity ratios of molecular lines are not uniform 
in the CMZ, suggesting spatial variation of gas prop- 
erties in the region. Such variations are most obvious 
in the HCN-to-CO ratio and in the ratio between ^^CO 
and other less abundant molecules (Fig. [7|). The HCN- 
to-CO intensity ratio is higher at the chain of molecular 
complexes, 0.06-0.08 at the CO (3-2) peaks, than in 
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other parts of the CMZ, where the mean ratio is 0.04 in 
the central 20''. To put it differently, 83% of HCN(4- 
3) flux detected in the central 20'' (340 pc, FWHM) 
is from 3" (50 pc) diameter regions at the five contin- 
uum peaks while the fraction is only 30-44% for the 
four CO lines. The HCN emission is clearly more con- 
centrated in the five molecular complexes. The ratio 
map Fig. [T^c) shows that the ratio has a peak value 
of about 0.12 near the millimeter continuum peak No. 
2. HCN emission generally traces dense molecular gas 
because of the large dipole moment of the molecule 
and hence its large critical density. Multi-transition 
analyses of HCN in the CMZ of NGC 253 suggested 
that HCN molecules are subthermally excited and that 
the most plausible density of the HCN - emitting ^as is 
^Ho 10^ cm~^ (iKnudsen et all 120071 : iPadione et~all 
119971 : iJackson et al.l 11995). Thus the concentration of 
HCN(4-3) emission in the molecular complexes strongly 
implies a large concentration of dense gas there, possibly 
in the form of many dense gas cores embedded in giant 
molecular cloud complexes. HCN enhancement at the 
largest cloud complexes may have some additional causes 
such as infrared pumping of HCN (Carroll & Goldsmith 
[l981), excitation by electron collisions (jPickinson et al.l 
IT977), and spatial variation of chemical abundance. An- 
other spatial variation is seen between ^^CO and less 
abundant molecules. Specifically, the galactic nucleus 
is the strongest peak in the integrated intensity maps 
of ^^CO (2-1) and (3-2) while the two adjacent molec- 
ular complexes are stronger in less abundant ^^CO, 
C^^O, and HCN (see Figs. [3] and [7|). It is interest- 
ing to note that gas around the continuum peak No. 2 
again has a peak in the CO isotopologue ratios. The 
off-nuclear peaks also dominated, and often produced 
double-peak appearance, in previous high-resolution (< 
5") observations of no n-CO molec ules (HCO+ and HCN, 
IKnudsen etH I ("2007') ; H2S,iMinh et al.l ([200 21); SiO and 
H^0+, iGlrcia-Bur Tllo et al.l (I2OOOI ): CS, iPeng et al.l 
(jl996)). This difference between the nuclear and off- 
nuclear gas peaks may be because the gas temperature 
is higher at the nucleus as we will see in §3.41 from the 
line-to-continuum ratio. 

3.3.5. Other Line Features 

Our line data also show features that were less clearly 
seen in previous works. First, the asymmetry of^^CO 
emission about the major axis, pointed out in IS06L is 
more evident in the new data. For example, in the 
^^C0(2-l) map in Fig. [3l the contours at around 2000 
K km s~^ are clearly asymmetric about the major axis 
that runs approximately through the continuum peaks 
marked with crosses. The north-western side of the nu- 
cleus, which is the side near to us, shows sharper out- 
ward decline in -'^^CO intensity than the south-eastern 
side , which is the far side of the galaxy. As inferred in 
lS06l . this is most likely due to high opacity and area fill- 
ing factor of ^^CO-emitting gas because the asymmetry 
is seen only in ^^CO. The asymmetry in optically thick 
lines can be caused by the vertical thickness of the CMZ 
and the plausible temperature gradient of molecular gas 
(i.e., lower temperature at the larger radii). It is be- 
cause, for a disk with thickness 2h and inclination i, the 
center of the disk surface is offset by hsini to the far 



side from the center of the disk volume when projected 
onto the sky. The former center is hottest in the pho- 
tosphere of an optically thick line while the latter is the 
true center of the disk and the hottest point for an op- 
tically thin emission. The offset can be enhanced if the 
molecular disk has the shape of a bic oncave lens, which 
is observed in the Galactic CMZ fe.g. JBallv et al.ll2010f ). 
The starburst wind from the CMZ of NGC 253 (e.g., 
iLehnert HeckmanI [T996h may be a reason for such a 
gas distribution. Similar near-side/far-side asymmetry in 
moleculai^ lines has been reported in the centers of NGC 
1068 (|Baker &: Sco ville 1998 |). IC 342 dMeie r & TurneJ 
I2M1I), NGC 1365 j[SakamotQ et al.ll20'07h . and NGC 4826 
(jBoone et al.ll2011l ): the last paper has more of possible 
cases and a model of the opacity effe ct. 

Second, one of the bubbles noted in lS06l . called SB2, is 
seen in our 1^00(2-1) map at (Aa, A^) - (20", 8") with 
respect to the radio source TH2 (see Fig. [3^), while SBl 
at around (—9", —7") is not evident in the integrated in- 
tensity map. The gas distribution at the position of SBl 
is indeed disrupted in both channel maps and the inte- 
grated map, but, unlike SB2, the feature does not appear 
to be semi-circular in our new data. The nature of SBl 
thus remains uncertain; it can be an expanding bubble 
in a complicated shape or other gas kinematical feature 
that looked like a bubble in lo w resolution. The third 
possible bubble inferred in lS06l from position- velocity di- 
agrams is identified with a depression at (9'.' 5, 3'.' 5) in the 
integrated intensity map (e.g.. Fig. [3fc). 

Third, the velocity map in Fig. [5] confirms the non- 
circular motion of molecular gas in the CMZ since the 
isovelocity contours run almost parallel to the major axis 
of the galaxy at P. A. ~50°. Our data also reveal sys- 
tematic deviation of the velocity field from the rest of 
the CMZ in the central few arcsec around the nucleus 
(or TH2). Velocity contours are kinked in that region 
in such a way that the kinematical major axis approxi- 
mately agrees with the galaxy major axis. This gaseous 
velocity field, in particular the strong noncircular mo- 
tion in the CMZ, is in sharp contrast to the circular 
velocity field of stars observed in the same region by 
iMiiller- Sanchez et al.l ([2010 ). Because of this large non- 
circular motion (i.e., lack of the usual spider pattern 
in velocity contours), the current data do not distin- 
guish whether the dynamical center of the galaxy is at 
the radio sour ce TH2 or ^ 0'.'7 from it as estimated by 
.Miiller- Sanchez et all (I2010 D. 

3.4. Comparison of Line and Continuum 

The molecular cloud complexes as well as their bright- 
ness temperature peaks are spatially associated with the 
'^1 mm continuum peaks and also with the centimeter- 
wave and infrared sources in and around them ( §3.2.3p . 
However, there are sometimes slight offsets of < 1" be- 
tween line and continuum peaks. This can be seen in the 
line emission maps in Figs. [3] and |4] where the millime- 
ter continuum peaks in Fig. [2] and Table [3] are plotted 
as crosses. For example, the 1 mm peak No. 2 is offset 
about O'.'5-l" west or northwest from the nearest peak 
in our line emission maps. This is seen even in optically 
thin C^^O as a 0'.'5 offset in both Fig. [3]i and Hi and is 
therefore unlikely to be due to line opacities. The line- 
continuum offset is less significant for the peak No. 3 
at the nucleus. A possible reason for the line-continuum 
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offset is that the dust emission fohows the heating from 
massive stars more truly than the molecular line emis- 
sion does because dust can absorb UV radiation and be 
heated while molecules are dissociated by the same ra- 
diation. If this is the case the offsets imply that active 
star formation occurred at the periphery of the molecular 
complexes or that the stars formed in groups are drifting 
away from the molecular complex. 

Regarding the intensity comparison between line and 
continuum emission, we have one note each for optically 
thick and thin lines. For optically thick ^^CO, the line- 
to-continuum ratio, or line equivalent width, is roughly 
uniform among the peaks No. 2 through 5 and is smaller 
than the mean value of the CMZ. (The peak No. 1 is 
too weak for this analysis.) The CO (3-2) equivalent 
widths are (3-4) x 10^ km s~^ while the mean value in 
the c entral 15'' is (15 ±3) x 10^ km s~^ (jSakamoto et al.l 
I2QQ8L and references therein). The smaller equivalent 
widths in these peaks are consistent with t heir high col- 
umn densities. The simple LTE model in Sakamoto et al. 
(2008) suggests column densities of TVhs/AV 10^^ 
cm~^(km s~^)~^, which agrees within a factor of a few 
with our estimate in §3.3.21 

For optically thin CO lines, the line-to-continuum ratio 
shows striking variation among the five peaks in contrast 
to the case of opaque ^^CO. The peak No. 3 at the nu- 
cleus is the brightest in the continuum but weaker than 
other peaks in ^^CO a nd C ^^O, as already noted in our 
previous obs ervat ions ()S06l ). One of the possible reasons 
discussed in lS06l was significant contribution from syn- 
chrotron or free-free emission to the continuum. This is 
not supported by the newly measured spectral index of 
the continuum at this peak, about 3, that is comparable 
to those at other peaks. A remaining plausible explana- 
tion for the stronger continuum at the nucleus is a higher 
temperature than in the other peaks. A higher tempera- 
ture reduces the equivalent widths of optically thin lines 
but does not change the equivalent width of a thick line 
in the LTE model (|Sakamoto et al]l2008[ ). The higher 
temperature at the nucleus is thus consistent with our 
observations of CO isotopologues. Presence of molecular 
gas at > 100 K in the central region o f NGC 253 has been 
inferred by observations of HC^ N (|Mauersberger et al.l 
p90), NHs (Mauersberger et al. 2003), high- J CO lines 
fBave t et al. 2004). and the combin ation of CO isotopo- 
logues and HCN (jMatsushita et aT]l2010[ ). A large frac- 
tion of such hot molecular gas may be at the nuclear 
peak (i.e.. No. 3) causing the low equivalent width of 
optically-thin CO lines there. Another remaining model 
is the selective photodissociation of ^^CO, C^^O, and 
other molecules showing less pronounced peaks at the 
nucleus (No. 3) than at No. 2 and 4. This may accom- 
pany the possible higher temperature and hence stronger 
heating radiation at No. 3. 

3.5. Comparison with the Center of our Galaxy 

It is interesting to compare the central molecular zone 
of NGC 253 with that of our Galaxy to learn the cause 
and consequence of starburst because the two are (among 
the) nearest starburst and non-starburst galactic cen- 
ters, respectively. Both have near edge-on viewing an- 
gles and owe their gas distributions to a galactic stellar 
bar (Scoville et al.. 1985; .Sorai et al.. .2000. : .Binnev et al.. 



119911 : iSawada et alll2004f ). The Galactic CMZ is the 
central half kpc region of high gas concentration with 
a total gas mass of /^5 x 10^ Mq (jP ierce-Price et al.l 
120001 : iDahmen et"allll998[ ) and a total infrared (A = 1 
-1000 /im) luminosity fro m interstellar dust of 4 x 10^ 
Lq (Sodroski et al.lll997[ ). Thus the starburst CMZ of 
NGC 253 has about an order of magnitude more gas and 
is about two orders of magnitude more luminous than 
the central region of our Galaxy (see references in 21 • 
A previous comparison of the two CMZs was made by 
Paglione et al. (1995) in IICN(l-O) at --50 pc resolution. 

We show CO (2-1) maps of the center of our Galaxy 
in Figure [8] for compa rison with NGC 253. The data 
were obtained by Saw ada et al.l (|2001[ ) using the Uni- 
versity of Tokyo 60 cm telescope and have the same 20 
pc resolution as our NGC 253 data. Two morphologi- 
cal similarities are evident from the comparison. First, 
both galaxies have a similar size of about 0.5 x 0.1 kpc^ 
at the half maximum intensity of the CO (2-1) emission. 
The two central molecular zones therefore have about the 
same size. Second, both galaxies show several prominent 
peaks of 20-50 pc sizes in their CMZs. The four most 
prominent CO peaks in the Galactic CMZ are Sgr A, B2, 
and C cloud complexes at / = 0°, 0.7°, and —0.5°, respec- 
tively, and the / = 1.3° cloud complex. About one third 
of the total C0(2-l) flux in the central 340 pc (FWHM) 
of our Galaxy is from the four 50 pc diameter regions 
centered at these molecular cloud complexes. They show 
peaks in the map of peak brightness temperature (Fig- 
ure [Sb). They also have associated peaks of centimeter 
and (sub) millimeter continuum emission as in the case of 
NGC 253. In particular, the Sgr B2 cloud complex con- 
tains the most prominent sites of star formati on in the 
Galactic center (Pi erce-Price et al.l 120001 : ISchuller et al.l 
2009; Bally et al. 20TqI and references therein). The four 
cloud complexes in the Galactic center have higher HCN- 
to-CO intensity ratios, or higher gas densities, than the 
area surrounding them according t o the J=l-0 line ob - 
servations of the Galactic center by Jackson et al.' (19 96'). 
This coincides with what we saw in NGC 253 ( ^3.3.4p and 
constitutes the third element of similarities between the 
Galactic center and the center of NGC 253, that is, the 
massive molecular cloud complexes in both CMZs have 
high HCN to CO ratios probably due to their high densi- 
ties. This agrees with what Paglione et al. (1995) found 
using HCN and CO (1-0) data of 3'' and 7'' resolution, 
respectively. 

The largest difference between the two galaxies is that 
the most prominent molecular cloud complexes in the 
center of NGC 253 have a factor of 3-4 higher C0(2- 
1) intensities, in both integrated and peak values, than 
those in the center of our Galaxy. The difference in the 
peak integrated intensities suggests, if the same CO-to- 
H2 conversion factor applies to both, that the cloud com- 
plexes in NGC 253 have higher masses and higher gas 
column densities at 20 pc scale than their Galactic coun- 
terparts. Since the molecular cloud complexes in the two 
CMZs are similar in size according to our CO (2-1) maps, 
the higher column density in NGC 253 would also mean 
higher mean volume density in the cloud complexes in the 
starburst CMZ. A caution here is that the uncertainty of 
the CO-to-H2 conversion factors in the centers of the two 
galaxies directly affects this comparison. We estimated 
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in ^3.3.2l tlie mass and mass surface densities of molecular 
gas in the center of NGC 253 using two methods inde- 
pendent of the conversion factor. Those estimates agree 
within a factor of a few to the mass derived from the 
adopted conversion factor of Xco = 1 x 10^^ cm~^ (K 
km s~^)~^, suggesting the mass uncertainty of 0.5 dex. 
The conversion factor in the cent er of our Galaxy also 
has its uncertainty. For example, iDahmen et al.l ()1998f ) 
mapped the Galactic center in optically thin C^^0(1- 
0), derived the mass of molecular gas, and compared the 
derived mass from other methods. Their weighted best 
estimate for molecular mass in the central 600 pc of our 
Galaxy is 4 times smaller than that from the CO (1-0) 
luminosity and the conve rsion factor we us e d for NGC 
253. On the other hand, iPierce-Price et al] (|2QQ0[ ) con- 
cluded from their Galactic center mapping in optically 
thin dust emission at 450 and 850 /im that the total gas 
mass in the cen tral 400 pc is 80% la rger than the best 
estimate mass of lDahmen et al.l (|1998f ) in the central 600 
pc. If this were the correct mass the conversion factor in 
the Galactic center is probably within a factor of 2 from 
what we adopted for NGC 253. Further evaluation of 
the conversion factors in two galaxies (and their possible 
variation within each galaxy) is beyond the scope of this 
paper, but readers are cautioned about this uncertainty 
in translating the observed difference in CO integrated 
intensities to difference in mass surface densities. 

The difference in the peak beam-averaged temperature 
is such that it is 50 K in NGC 253 and 17 K in the 
Galactic center at the same 20 pc resolution in the same 
CO (2-1) line. The peak brightness temperature in the 
Galactic CMZ is lower even at higher spatial resolution 
and with less beam dilution. The peak Tb of CO (1-0) 
and (3-2) emission in the Galactic CMZ is a lmost always 
below 30 K at 0.8 pc resolution (Oka et al.l [19981 12007). 
Only in the Sgr A and B regions lOka et al.l (I2001D found 
a cloud whose peak temperature exceeds 30 K; the peak 
Tb is 38 K and 31 K, respectively, at 2.5 pc resolution. 
We therefore suggest that the difference in peak Tb at 
20 pc scale is mainly because the molecular complexes 
are warmer in NGC 253 than in the Galactic CMZ, al- 
though 50 K gas is not excluded in the Galactic CMZ at 
sub-parsec scales or behind the CO photospheres. The 
warmer gas in NGC 253 is consistent with the larger lu- 
minosity of the starburst nucleus and in accordance with 
the higher dust temperature there. The dust tempera- 
ture in the Galactic CMZ has been estimated to be ~20 
K from SEP fitting in the subm illimeter-to-far infrared 
bands ([Fierce-Price et all I2000L and references therein) 
while dust in the center of NGC 253 wa s estimated also 
from SEP fitt ing to be a t or above 3 K ("Radovich et alj 
HOOII iMelo et al. 2001 iWeifi et al.l 12008). The high- 
temperature gas in the center of NGC 253 has also been 
pointed out from the line excitation analyses mentioned 
in ^3.41 Our Tb observations add a new temperature 
constraint that does not involve such delicate issues in 
the line excitation modeling as the assumption on cloud 
velocity structure. Finally, the higher density and tem- 
perature of the molecular gas in the center of NGC 253 
suggest that gas thermal pressure is also highe r there 
than in the Galactic CMZ, where ISawada et all (|200lD 
estimated P/k = nu^T ~ 10^ cm~^ K. The difference 
will be about an order of magnitude on average if the 



density increase in NGC 253 is in the microscopic nn^. 
4. SUMMARY 

The center of the nearby starburst galaxy NGC 253 
has been imaged in five molecular lines and continuum 
at 0.87 mm and 1.3 mm at resolutions of 20-40 pc. A 
comparison was made in CO (2-1) with the center of 
our Galaxy at the same linear resolution. Our main ob- 
servations are the following. 

1. Molecular gas in the center of NGC 253 is con- 
centrated in an elongated region on the sky, or a highly 
inclined disk, with the total extent of ~1 kpc, as has 
been seen in lower-resolution observations. This region 
has a FWHM size of about 300 pc x 100 pc when ob- 
served in CO (2-1) at 20 pc resolution. This size is 
approximately the same as the FWHM size of the Cen- 
tral Molecular Zone (CMZ) of our Galaxy observed in 
the same line at the same linear resolution. 

2. Within the CMZ of NGC 253, there is a string of 
peaks in molecular lines and 870-1300 /im continuum 
emission. We reported the parameters of five most con- 
spicuous peaks of < (20-40) pc extent. They are most 
likely molecular cloud complexes comparable in size with 
such molecular cloud complexes as Sgr A, B2, C and the 
/ = 1.3° complex in the CMZ of the Milky Way. 

3. The molecular complexes in NGC 253 are warm, 
dense, and have high column densities. They have ^^CO 
peak brightness temperatures of ~50 K, molecular gas 
column density on the order of 10^ Mq pc~^, gas mass 
on the order of 10^ Mq within a size scale of 20 pc, 
and volume- averaged gas density of riHs ^ 4000 cm~^. 
HCN-to-CO ratio is higher at the molecular complexes 
than in the surrounding area, as was observed in the 
Galactic CMZ, plausibly suggesting higher gas densities 
in the cloud complexes. There are indications from line 
and continuum ratios that the molecular complex at the 
dynamical center has different properties, likely a higher 
gas temperature, compared to other clouds. Compared 
with their counterparts in the center of our Galaxy, the 
molecular complexes in the CMZ of NGC 253 have peak 
brightness temperatures and column densities that are 
about a factor of 3 higher at the same spatial resolution, 
although the comparison of column density is subject to 
the uncertainty in the conversion factor. 

4. The molecular complexes are spatially associated 
wi th the centime t er- wa ve radio peaks such as those found 
bv lTurner fc Hoi (|1985f ). This coincidence, as well as the 
coincidence of the molecular complexes and other star 
formation tracers, suggests that these molecular com- 
plexes are the sites of the most intense star formation 
in the center of NGC 253. This is in parallel to the Sgr 
B2 cloud complex being the most prominent site of star 
formation in the Galactic CMZ. 

5. We find signs that the high opacity of ^^CO pre- 
vents it from truly tracing gas distribution both in indi- 
vidual cloud complexes (because of self-absorption and 
cloud overlapping) and in the overall structure of the 
CMZ (because we only see the warm surface of the CO 
photosphere.) Caution is needed about this bias when 
analyzing the gas parameters using only ^^CO data. 
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Table 1 

Log of SMA observations 



observation date 


main line 


/lo 


-^track 


A^ant 


T225 


{rsys,DSB> 


-^baseline 


Tint 






GHz 








K 


m 


hr 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


2004 Aug.-2005 Nov. 


CO(2-l) 


224.650 


6 


7 or 8 


0.07-0.12 


130 


8-509 


5.7 


2004 Sep. 


CO(3-2) 


340.508 


1 


7 


0.04 


279 


7-124 


4.2 


2004 Nov. 


HCN(4-3) 


349.222 


1 


7 


0.09 


447 


7-124 


3.3 



Note. — (3) Frequency of the local oscillator. The center frequencies of upper and lower sidebands 
are /lo i 5.0 GHz. Each sideband is 2 GHz wide. (4) Number of tracks (i.e., nights) used. The 1.3 mm 
observations used array configurations called compact, compact-north, extended, and very-extended for 1, 
2, 2, and 1 night (s), respectively. (5) Number of participating antennas. (6) Zenith opacity at 225 GHz 
measured at the Caltech Submillimeter Observatory next to the SMA. Opacity in the 350 GHz band is 
approximately 3.5r225- (7) Median of double-side-band system temperature toward NGC 253. (8) Range 
of projected baseline length for the galaxy. (9) Total integration time on our target positions. The short 
Tint of our 230 GHz observations for the number of tracks is because some off-center positions were also 
observed in these tracks. 



Table 2 

Data parameters 



line, band 


resolution 




rms noise 




/sD 


^SD 


SD ref. 




" X " 


km s"-*^ 


mJy beam~^ 


K 


% 


II 




12C0(2-1) 


1.14 X 1.13 


10 


39 


0.69 


79 


23 


1 


13C0(2-1) 


1.62 X 1.49 


20 


17 


0.18 


89 


23 


1, 2 


Ci80(2-1) 


1.62 X 1.49 


20 


18 


0.19 


97 


23 


1 


1.33 mm continuum 


1.10 X 1.10 




2.5 


0.050 


>64 


~30 


3, 4 


i2CO(3-2) 


1.82 X 1.10 


10 


215 


1.2 


62 


22 


5 


HON (4-3) 


2.20 X 1.71 


30 


130 


0.34 


72 


20 


6 


0.87 mm continuum 


1.80 X 1.10 




23 


0.12 


60 


15 


7, 8 



References. — 1. IHarrison et alj ((19991): 2. Ilsrael Baai (1200^): 3. IKriigel et al.l 
(1990); 4. Mauersberger et al. (1996); 5. IBavet et al.l (|2004l ): 6. IXnudsen et al.l (|2007l ): 7. 
Alton et al. (1999); 8. Seaquist et al. (200l) 

Note. — The resolutions are in full width at half maximum (FWHM). \" = 17 pc. 
The last three columns are fractions of single dish (SD) fluxes detected in our data, the 
FWHM beam sizes of the single dish observations, and the sources of the single dish data. 
The comparisons are made at the center positions given in the single dish papers. Line 
contamination in the single-dish continuum observa tions has b een subtracted adopting the 
estimates of -"^^CO contamination by Mauersberger et al.l (|199 6) and Seaquist et al. (2004); 
specifically the subtraction was 25% and 30% for the 1.33 and 0.87 mm, respectively. The 
uncertainty of /sd is about 15% and 20% for 1.3 mm and 0.8 mm data, respectively. 
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Table 3 

Near 1 mm Continuum Peaks in the Center of NGC 253 



No. 


OL 


b 


/l.3 


^0.87 


^0.87/1.3 




^9 




association 




(00^47™) 


(-25°170 


Jy beam ^ 


10^ 


Mo pc-2 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) (10) 


(11) 


1 


33! 63 


13(^2 


0.017 


(0.04) 




0.9 


(0.2) 


0.4 (0.2) 


H20_2 


2 


33? 30 


15^6 


0.064 


(0.22) 


2.7 


3.5 


(1.3) 


1.5 (0.9) 


THl 


3 


33! 17 


\r!\ 


0.082 


0.33 


3.0 


4.5 


2.0 


1.9 1.3 


TH2?^nucleus, TH4=H20_1, TH3,5,6 


4 


32! 99 


19(^7 


0.063 


0.27 


3.6 


3.4 


1.6 


1.4 1.0 


TH7(=SSC),TH8 


5 


32! 84 


21(^3 


0.044 


0.21 


3.5 


2.4 


1.3 


1.0 0.8 


TH9 



Note. — Columns 2 and 3 : Peak position of 1.3 mm continuum in J2000. Seconds of right ascension and arcseconds of 
dechnation, respectively. Columns 4 and 5 : Peak flux density of 1.3 mm and 0.87 mm continuum, respectively. The conversion 
factor to brightness temperature is 19.9 K/(Jy beam"-*^) for 1.3 mm and 5.2 K/(Jy beam"-*^) for 0.87 mm. The 0.87 mm intensities 
in parenthesis are measured at the corresponding 1.3 mm peak positions because 0.87 mm data lack an isolated peak. Column 
6 : Ratio of 0.87 mm to 1.3 mm continuum flux densities measured in the data with matching u-v coverage and 2^^ resolution. 
The fractional error of the ratio is 20%. The peak No. 1 was too weak at 0.87 mm to reliably derive the ratio. Columns 7 and 
8 : Peak gas surface density estimated from 1.3 mm and 0.87 mm data, respectively, with = 5.4 x 10^ x (/1.3/Jy beam~^) 
Mq pc~'^ and = 6.0 x 10^ x (/0.87/Jy beam"-*^) Mq pc~^. See ^3.2.2l for assumptions and uncertainties. Columns 9 and 10 : 
Gas mass of each peak in the 1.3 mm and 0.87 mm beam, respectively. The beam FWHM is 19 x 19 pc^ and 31 x 19 pc^ in our 
1.3 mm and 0.87 mm data, respectively. Column 11 : Association with sources known in other wavelengths. Compact centimeter 
sources of iTurner H o (1985) are called THn, among which TH2 is thought to be at the galactic center (or within 2^^of it; 
IFernandez-Ontiveros et al. 2009). The nuclear region has two water masers H20_l (at Vlsr ~ 120 km s"-*^) and H20_2 (170 
km s~^) associated with star formation (Henkel et al. 2004; H ofner et al.ll2006l : IBrunthaler et al.i i2009). H20_2 is likely associated 
with the molecular cloud complex that has the peak line velocity of 170 km s~^ (see Fig. [6]) and includes the continuum peak No. 
1 listed above. SSC is the most luminous mid-IR source in the nucleus and is suggested to be a super star cluster ()Keto et al.l 
[19991 



Table 4 

Molecular-cloud Complexes in the Center of NGC 253 



No. 

(1) 


00^47™ 
a 
(2) 


-25°17' 
6 

(3) 


lAV 

(4) 


CO(2-l) 
Tp I'AV 
(5) (6) 


CO(3-2) 
lAV Tp 
(7) (8) 


13 

lAV 

(9) 


CO(2-l) 
(10) (11) 


C 
lAV 
(12) 


180(2-1) 
(13) (14) 


HCN(4-3) 
lAV Tp 
(15) (16) 


1 


33! 67 


13(^1 


5.5e3 


52 


4.9e3 


4.2e3 


41 


5.1e2 


6.0 


0.11 


1.4e2 


2.2 


0.028 


2.6e2 2.9 


2 


33! 38 


15(^8 


6.1e3 


43 


5.3e3 


4.5e3 


36 


8.5e2 


10.4 


0.17 


2.8e2 


4.7 


0.054 


3.7e2 3.6 


3 


33! 18 


l/.'4 


6.5e3 


34 


5.6e3 


5.7e3 


34 


7.3e2 


5.1 


0.14 


2.2e2 


1.5 


0.040 


4.6e2 2.5 


4 


32! 99 


19(^6 


6.2e3 


48 


5.6e3 


5.2e3 


40 


8.4e2 


7.7 


0.16 


2.3e2 


2.8 


0.042 


4.1e2 5.1 


5 


32!82 


20'.' 9 


5.0e3 


47 


4.6e3 


3.7e3 


36 


7.9e2 


10.5 


0.19 


2.1e2 


4.5 


0.047 


2.6e2 3.0 



Note. — Columns 2 and 3 : Peak positions of -"^^ CO (2-1) integrated intensity, except for peak No. 5 whose position is 
from the -'^'^ CO (2-1) integrated map because -'^^00(2-1) does not show a peak. Seconds of right ascension and arcseconds 
of declination, respectively, in J2000. Columns 4, 7, 9, 12, 15 : lAV is the integrated intensity in K km s"-*^ at the position 
of (2) and (3). Columns 5, 8, 10, 13, 16 : Tp is the peak intensity in K at the position of (2) and (3). Both lAV and Tp 
are measured from the data at their native spatial resolutions in Table O Column 6 : I' AV is the integrated intensity in 
K km s"-*^ at the position of (2) and (3) measured after the data are convolved to the spatial resolution of the ■'^^00(2-1) 
and C-'^^0(2-l) data. Columns 11 and 14 : Mean opacity of the line, f, estimated from the ratio of integrated intensities 
between the isotopologue and ^-^CO. 
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4 

g 3 

1 




: co(2-i) 




^ 1 .86+04 Jy km/s [ 1, 




1 .Oe+03 K km/s^^Hl 







100 200 300 400 
Velocity [km/s] 



CO(3-2) 
2.3e+04 Jy km/s 
5.8e+02 K km/s 




100 200 300 400 
Velocity [km/s] 




200 300 400 
Velocity [km/s] 



HCN(4-3) 
0.15- 9.6e+02 
2.3e+01 




200 300 400 
Velocity [km/s] 




HNCO(10o,o-9o,9) 




100 200 300 400 
Velocity [km/s] 



100 200 300 400 
Velocity [km/s] 



Figure 1. Line emission at the center of NGC 253. The spectra are taken from the SMA data cubes after correcting the cubes for the 
primary beam attenuation and convolving them to 20'' resolution (FWHM). The measured line fluxes are given in the legend of each panel 
in units of Jy km s~^ (second line) and K km s~^ (third line). The continuum flux density at the center at the same resolution is 0.61 Jy 
beam~^ (= 35 mK) and 1.2 Jy beam~^ (= 30 mK) respectively at 1.3 mm and 0.87 mm. No correction for missing flux is made. 
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20 40 60 80 




100 200 300 




10 5 -5 -10 

R.A. offset [arcsec] 

Figure 2. Radio continuum in the center of NGC 253. (a) 1.3 mm, (b) 0.87 mm, and (c) 1.3 cm and other sources for comparison. In (a) 
and (b), contour steps are 7.5 mJy beam"-*^ (= 0.05 K = 3cr) and 40 mjy beam"-*^ (= 0.21 K = 1-7(7), respectively, and peak intensities are 
82 mJy beam"-*^ (= 1.6 K) and 0.32 Jy beam"-*^ (= 1.7 K), respectively. The five brightest peaks of 1.3 mm continuum are marked with 
plus signs in (a) and (b) and with circles in (c). Panel (c) shows the 1.3 cm compact sources listed in Ulvestad &; Antonucci (1997, Table 
6) as diamond symbols, each of which has an area proportional to its fiux density. The brighter sources found earlier by Turner Hd 
{1985) are indicated by red squares with THn names while others are shown in magenta. TH2 is the brightest source at 1.3 cm. The green 
asterisk is the most IR-luminous star cluster in this region and is associated with TH7 (Keto et al. 1999; Fernandez-Ontiveros et al. 200^. 
Two blue plus signs are H2O masers, one of which is at TH4 (Henkel et al. 2004; Hofner et al. 2006). The origin of the offset coordinate 
in this and subsequent figures is q;=00^47™33!182, (5=-25° 17^17(^148 (J2000), the position of TH2 measured bv .Lenc &: Tingavi (|2006i ). 
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Figure 3. Molecular lines in the center of NGC 253. The line name and the synthesized beam (FWHM) are shown in the 
top-right and bottom- left corners, respectively. The five plus signs are at the positions of the 1.3 mm continuum peaks. The 
dotted lines are the half-power contours of the (mosaicked) primary beams, for which the maps are corrected. The dashed lines 
mark the area shown in Figs. [2]and|4l The unit of integrated intensity is K km s~^. The n-th contour is at 46 x n^, 78 x n^"^, 
17 X r^^•^ 18 x •^ and 36 X n ■ , respectively, from (a) to (e). 
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10 20 30 40 




2 4 6 8 10 12 
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10 5 -5 -10 

R.A. offset [arcsec] 

Figure 4. Peak brightness temperature maps of molecular lines in the center of NGC 253. The Ray leigh- Jeans brightness temperatures 
of the lines (in excess of any continuum) are in units of K. The five plus signs are at the positions of the 1.3 mm continuum peaks. 
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20 10 -10 -20 
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Figure 5. Mean velocity map of ■'^^00(2-1) emission in the center of NGC 253. Contours are in 25 km s ^ steps. The five plus sign s are 
at th e positions of the 1.3 mm continuum peaks. The approaching line of nodes of the galaxy is at the position angle of 51° (^Jarrett et al.l 
[2001) . 
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Figure 6. Line profiles at the five most prominent molecular cloud complexes in the galactic center. Each spectrum is from the cloud 
position in Table 2] and is corrected for the primary beam attenuation. The variation of line shapes between different molecular lines, seen 
in the molecular complexes No. 2, 3, 4, and 5, suggests self-absorption of -"^^CO or multiple velocity components. 
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Figure 7. Ratio of line integrated intensities, (a) i3cO(2-l) to i2cO(2-l). (b) 01^0(2-1) to i2cO(2-l). (c) HCN(4-3) to CO(3-2). For 
each ratio, the two integrated intensity maps were convolved to the same resolution and expressed in units of K km s"-*^ before calculating 
the ratio at locations where both lines are detected above 5 a. Plus signs show the locations of the five 1.3 mm continuum peaks. 
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Figure 8. The cen ter of the Milky Way galaxy in -'^^00(2-1) for comparison with the center of NGC 253. The single-dish data are from 
ISawada et al.l (|2001l) and have a resolution of 23 pc (FWHM) at the Galactic center (distance 8.5 kpc). The dashed lines mark the same 
linear-scale area as they do in Fig. [3] (a) Integrated intensity. The contours are at 46 x n^-^ K(Tmb) km s"-*^. The brightest peaks at / = 
0°, 0.7°, -0.5°, and +1.3° are the Sgr A, B2, C clouds and the I = 1.3° cloud, respectively, (b) Peak brightness temperature with 2 K 
contour intervals. Peak intensities were searched only at |Vlsr| > 30 km s"-*^ to exclude emission from the disk and arms in the foreground 
or background of the Galactic center. 



